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A representation-theoretic approach to special functions was developed in the 
40-s and 50-s in the works of I.M.Gelfand, M.A.Naimark, N.Ya.Vilenkin, and their 
collaborators (see [V],[VK]). The essence of this approach is the fact that most clas- 
sical special functions can be obtained as suitable specializations of matrix elements 
Q\ ' or characters of representations of groups. Another rich source of special functions 

is the theory of Clebsch-Gordan coefficients which describes the geometric juxtapo- 
sition of irreducible components inside the tensor product of two representations (cf. 
[VK]). Finally, in recent works on representations of (quantum) affine Lie algebras 
Q ' it was shown that matrix elements of intertwining operators between certain repre- 

, sentations of these algebras are interesting special functions - (q-)hypergeometric 

7-H ' functions and their generalizations [TK, FR]. 

In this paper we suggest a general method of getting special functions from 
representation theory which unifies the three methods mentioned above and allows 
one to define and study many new special functions. We illustrate this method by 
' a number of examples - Macdonald's polynomials, eigenfunctions of the Sutherland 

7-H , operator, Lame functions. Other examples will be described in our future papers. 

1. Vector-valued spherical functions. Let TL be a Hopf algebra over C, 
• H C TL be a subgroup of group-like elements, V, W, U be irreducible H-modules. 

\ Let $ : V — > W ® U be an intertwining operator for TL. From this data one can 

construct two kinds of special functions. 

1. Vector-valued matrix elements. Set f vw $(h) —< w,<&hv >E U , h E H , 
(-} • v E V,w E W* . We call this function on H a [/-valued matrix element. 

2. Vector-valued characters. Let V — W. Set x$(/i) = Tr |y ($/i) eU, h E H. 
' We call this function on if a [/-valued character. 

Example 1. If U = C, V = W, $ = Id, and we get the classical matrix 
k>< | elements and characters. 

Example 2. The numbers < f vw <s>(e),u >, u £ U* , where e is the identity 
element in H, are the Clebsch-Gordan coefficients for Ti (we assume that v, w, u 
are taken from some bases of V, W*, U*). 

Example 3. Let TL be the quantum affine algebra U g (g) corresponding to a 
finite -dimensional simple Lie algebra g (q can be equal to 1), H = {e}. Let V, W 
be Verma modules over TL with the same central charge, and let U = Ui(zi) ® • • • ® 
Un{z n ), where Zi E C*, Ui are finite-dimensional representations of 7i, and Ui(z) 
denotes the representation of TL in the space Ui defined by nudz) ( a ) = z Acg ^-Ku i (a), 
a E TL, deg(a) being the homogeneous degree of a. Let $ : V — > W ® U be 
a vertex (=intertwining) operator (here (g> should be understood as a completed 
tensor product). Let v, W be the vacuum vectors of V,W*. (For more details 
see [FR]). Then the Ui <g> • • • <g> C7 n -valued matrix element f vw &(e) is a function of 
n complex variables Zx,...,z n . For q — 1 such functions arise in conformal field 
theory as correlation functions on the sphere; for q ^ 1, they appear as correlation 
functions of solvable lattice models in statistical mechanics. It is known that they 
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satisfy the (classical or quantum) Knizhnik-Zamolodchikov (KZ) equations and 
express via hypergeometric (q-hypergeometric) functions and their generalizations. 
Thus, solutions of (quantum) KZ equations can be regarded as Clebsch-Gordan 
coefficients for (quantum) affinc Lie algebras. 

2. A generalized Peter- Weyl theorem. In this section H = G is a finite 
group or a compact Lie group, and H is the convolution algebra of (generalized) 
functions on G. Notation: dg is the invariant probability measure on G; R(G) is 
the set of irreducible unitary representations of G; U G R{G); (-,-)u is the inner 
product in U; L 2 {G,U) is the Hilbert space of U- valued L 2 -functions on G with 
respect to dg, with the inner product < /i,/2 >= J G (fi(g), f2(g))udg; L 2 (G,U) G 
is the subspace of L 2 (G, U) spanned by all functions \ such that xihgh^ 1 ) = h\{g), 
g, he G. 

For every V G R{G), fix an orthonormal basis By of V. For any V, W G 
R(G), let Xyw be the space of G-homomorphisms: $ : V — > W ® t/. This is a 
Hilbert space since it is isomorphic to the space of G-invariants in V* <g> W ® [/. 
Let ByvF be an orthonormal basis in Xyw- Let G L 2 (G,U) be defined by 

fvw$(g) = (w,<&gv)w E U, g e G,v e V,w E W,Q e X vw . We have f V w<s>(gh) = 
f(hv)w<b{g) = gfv( g -iw)$(h)- Also, define the functions x<s>(g) = Tr |y (<$>g) = 
Y,veBv f™*{g), $ G ^vy- Clearly, X * G L 2 (G, C/) G . 

Theorem. fij {(dimT/dimVF) 1 / 2 /^^* : w G G G B VW ,V, W G 

i?(G)} is an orthonormal basis of L 2 (G, U). 

(ii) {x$ : $ G V G R(G)} is an orthonormal basis of L 2 (G, U) G . 

When U — C, this theorem reduces to the standard Peter- Weyl theorem [K]. 

Remark. This theorem can be generalized to the case of quantum groups (cf. 
[EK1]). 

3. Macdonald's polynomials as vector-valued characters. Let G = 

SU(n), TL be the algebra of functions on G, and let H be the Cartan subgroup 
of G (diagonal matrices). Macdonald's polynomials P\(q,k,h) (k G Z, h G H, A is 
a dominant integral weight for G) are uniquely defined by the following properties: 

1. P\(q, k, •) are trigonometric polynomials on H symmetric under the action of 
the Weyl group S n . 

2. For fixed q and k, P\(q, k, •) are orthogonal on H with respect to the weight 
|A| 2 , where A(q,k,h) = n^ioIL^ 1 - q m e <a ' i> ), where £ G h = Licff is such 
that h = e^, and a runs over all positive roots of G. 

3. P\ = x\ + J2v<\ c >-vXvi where v is a dominant integral weight of G, and \v 
is the character of the irreducible representation of G with highest weight v, and 
c\ u are constants depending on q, k. 

Example: P x (q,l) = X\- 

In the special case q = 1, Macdonald's polynomials are called Jack's symmetric 
functions. 

Let V — L\ be the finite-dimensional irreducible representation of G with highest 
weight A. Let U = S kn C n , where A; is a positive integer, be a representation of G. 
Let $ be a nonzero intertwining operator V — > V^®^. Such an operator exists iff A > 
kp, where p is the half-sum of positive roots of G, and if it exists, it is unique up to a 
factor. Let W(h) be the Weyl denominator: W(h) = Ila>o( e ^ <a ' 5> _ e"^ <a,5> ), 
h = e^. Let A = v + kp, where v is any dominant integral weight. Define the 
functions tp v (k,h) — W(h)~ k Tr\ v These functions take values in the zero 
weight component of U, which is one-dimensional. Thus, we can regard them as 
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scalar functions, choosing the normalization in such a way that the coefficient to 
e <»,£> i n ^ i s i. 

Theorem. [EK1] The functions il> v (k,h) are the Jack's symmetric functions 
P v (l,k + l,h). 

Remark. The orthogonality of {Vv} immediately follows from the generalized 
Peter- Weyl theorem. 

In the special case G = SU(2) the polynomials P v (l,k + l,h) arc the classi- 
cal Gegenbauer polynomials - the even trigonometric polynomials in one variable 
orthogonal with respect to the weight sin 2/c+2 x. 

Let now H be the quantum group U q (sl n ). The rest of notation is as above 
(the modules L\ and S kn C n are the 7Y-modules obtained by q-dcformation of 
the corresponding 5J7(n)-modules). Let W Qtm (h) = Yl'^2oYla>o( < l me ^ <a ^ > ~ 
? -m e -i<a,i>) ; h = e | Dennc tne functions ip„(q,k,h) = W qik (h)- lr Tr\ v (<S>h), 
with the normalization described above. 

Theorem. [EK1] The functions ip„(k, h) are the Macdonald's polynomials P v {q 2 , k+ 
l,h). 

4. Diagonalization of the Sutherland operator. The Sutherland op- 
erator is the Hamiltonian of a quantum many-body problem (sec [OP]): H = 

+ CSi<j smn ( a; i — x j)~ 2 (C i s a constant). It is known [OP] that 
H has n functionally independent commuting quantum integrals L\, L2, L n : 
L m = YTj=i JE:+ lower order tcrms ( L i = YTj=i afj) L 2 = -2H), and [L i} H] = 
for all i. Therefore, it makes sense to consider the system of differential equations 
Liip = Aitp, Aj e C (the eigenvalue problem). For each set of eigenvalues Aj this 
system has nl linearly independent solutions. For special values of {Ai} one of 
these solutions expresses via Jack's polynomials, but in general solutions are tran- 
scendental and express via generalized hypergeometric functions. Here we interpret 
these solutions as vector- valued characters for gl n . 

Let H = U(gl n ), f) = C™ be the Cartan subalgebra in H = cxp(()). Let 

V be the Verma module over gl n with highest weight A G C". Let /1 e C be 
such that C = u,(u, + 1), and let W be the module spanned by all functions of the 
form z^ +mi 4 +m2 . . . z£ +m ", rrij £ Z, £\ mj = 0, with the action of g[„ given by 
Eij 1 — ► z i~§^- ~ M<% ((Eij)ki = S ik Sji). Let $ : V — > V ® W be an intertwining 
operator for gl n (<g> is the completed tensor product). <& is unique up to a factor if 

V is irreducible (which happens for a generic A). 

Define the functions (f>\(xi, x n ) = VF(e^) _1 Tr|y($e ? ), where £ = 2Y^j=i x jEjj 
(W is the Weyl denominator). As in Section 3, we regard <f>\ as scalar functions. 

Theorem. [E] (i) The function satisfies the system of equations Liipx — 
Aitp\, where Ai = pi(X + p), and p m are symmetric polynomials with the highest 
term proportional to ^ . A™ 1 . 

(ii) For a generic A, the functions ip a (\+p)- P : ° € S n , form a basis in the space 
of solutions of the system = pi(X + p)tp. 

Remark. This theorem can be generalized to the case of quantum Ql n , which 
provides a diagonal basis of functions for Macdonald's difference operators (cf. [Ch]) 

5. Lame functions as vector- valued characters of s^ 2 - We preserve the 
notations of Example 3 from Section 1. Let jj = SI2, V = W = M\, k is the Verma 
module with highest weight A and central charge k, n — 1, z\ — 1, and U\ is the 
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same as in Section 4. Consider the function 

F(x, t) = (Tr| Mfe/2 . fc ( e -^ d +^+iW)^ ~* TT\ v {<P e -^ 2Td+ ^+^), 

where d is the homogeneous gradation operator and h — diag(l, —1) est This is 
a 1-point correlation function of conformal field theory on the torus. As in Sections 
3,4, we regard F as a scalar function. 

Theorem. [EK] F satisBes the Schrddinger equation —2iri(k + 2)^ + = 
fi(fx + l)(p{x + J, t) + c)F, where p is the Weierstrass elliptic p-function, and c is 
a constant. 

If fi is a positive integer, it is possible to express F as a /x-dimensional integral of 
products of powers of theta- functions and exponents. When k — > — 2 (critical level), 
the above equation becomes Lame equation [WW] ; finding the asymptotics of the 
integrals, we recover classical formulas from [WW], expressing Lame functions via 
theta- functions (see [EK]). 
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